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Abstract 

This note provides explicit techniques to compute the exponentials of a variety of structured 
4x4 matrices. The procedures are fully algorithmic and can be used to find the desired 
exponentials in closed form. With one exception, they require no spectral information about 
the matrix being exponentiated. They rely on a mixture of Lie theory and one particular 
Clifford Algebra isomorphism. These can be extended, in some cases, to higher dimensions 
when combined with techniques such as Given rotations. 
PACS Numbers: 03.65. Fd, 02.10. Yn, 02.10. Hh 



1 Introduction 

Finding matrix exponentials is arguably one of the most important goals of mathematical 
physics. In full generality, this is a thankless task, [1]. However, for matrices with structure, 
finding exponentials ought to be more tractable. In this note, confirmation of this phenomenon 
is given for a large class of 4 x 4 matrices with structure. These include skew-Hamiltonian, 
perskewsymmetric, bisymmetric (i.e., simultaneously symmetric and persymmetric -e.g., sym- 
metric, Toeplitz), symmetric and Hamiltonian etc., Some of the techniques presented extend 
almost verbatim to some families of complex matrices (see Remark (3.4), for instance). Since 
such matrices arise in a variety of physical applications in both classical and quantum physics, 
it is interesting that their exponentials can be calculated algorithmically (these lead to closed 
form formulae), for the most part, without any auxiliary information about their spectrum. 
For general symmetric matrices, however, the spectral decomposition of a 3 x 3 matrix is 
needed [see, however, iii) of Remark (2.1)]. On the other hand, this spectral decomposition 
can itself be produced in closed form. Thus, even for such matrices the techniques described 
here can be justifiably called closed form methods. For brevity, this note only records explicit 
algorithms for finding these matrix exponentials - the resultant final formulae can easily be 
written down once the reported procedures are implemented. 

The methods discussed below are of two types. The first, which is more versatile, relies 
on an algebra isomorphism of real 4x4 matrices with H <g) H . This algebra isomorphism, 
known from the theory of Clifford algebras and which ought to be widely advertised, was 
used in a series of interesting articles by Mackey et al., [2, 3, 4] for finding eigenvalues of some 
of the structured matrices discussed here. The present note can be seen as a contribution 
of a similar type. It is emphasized that for the preponderance of the matrices, considered 
here, this algebra isomorphism alone is needed. In particular, in this note no use is made of 
any of the structure preserving rotations used in [2, 3, 4] ever - see ii) of Remark (3.2). The 
second is based on the observation that several "covering" space Lie group homomorphisms, 
when made explicit, contain in them a recipe for finding exponentials of matrices belonging 
to certain Lie algebras. This circumstance renders the exponentiating of some 4x4 matrices 
(real/complex) equivalent to the job of finding the exponential of 2 x 2 matrices - which can be 



done in closed form. This method is, however, applicable only to a limited family of matrices. 
Therefore, this method is presented in an appendix. 

It is worth noting that, though most of the structured matrices considered here were chosen 
for their importance in applications, the real enabling structure is that present in H ® H. 
This is especially illustrated by certain normal matrices [see Definition (3.1)]. 

The balance of this note is organized as follows. In the next section some notation and one 
observation which is used throughout is recorded. In the same section the relation between 
H ® H and gl(A,R) is presented. The third section discusses a wide family of matrices 
which can be exponentiated using the aformentioned algebra isomorphism. The final section 
offers conclusions. In the first appendix the second approach to exponentiating matrices in 
p(4, R) and so(2, 2, R) (see next section for notation) is presented in a manner that makes 
the connection to the covering space homomorphism between SU(2) x SU(2) and 50(4, R) 
explicit (see Remark (5.1)). In the final section thirteen classes of matrices are listed which 
can be exponentiated by mimicking verbatim two situations studied earlier. 

In closing this introductory section it is noted that by combining these techniques with 
techniques such as structure preserving similarites, e.g, Givens rotations, [5], one can extend 
these results, in many cases, to find algorithmically the exponentials of structured matrices 
of size bigger than four. In other words, one can use such similarities (normally used in the 
literature for reduction to canonical forms) to reduce the exponential calculation to dimension 
four or lower of matrices with similar structure. In principle, this would provide closed form 
formulae for the exponentials of such structured matrices, since one can explicitly write down 
the desired Givens type similarities. However, it is more accurate to say that this implies 
an algorithmic procedure for exponentiating such matrices. For matrices for which this is 
possible (e.g., symmetric matrices), the details of this procedure is routine and hence will not 
be pursued here. 



2 Notation and Preliminary Observations 

The following definitions and notations will be frequently met in this work: 
• gl(n,R) and gl(n,C) represent the real (resp. complex) n x n matrices. 



• sl(n,R) and sl(n,C) represent the real (resp. complex) traceless matrices. SL(n,R) 
and SL(n, C) represent the real (resp. complex) matrices of determinant one. 

• SU(n) represents the n x n unitary matrices of determinant one. su(n) represents the 
n x n skew-Hcrmitian, traceless matrices. 

• R n represents the matrix with 1 on the anti-diagonal and zeroes elsewhere. p(n, R) and 
p(n, C) represent the nxn real (resp. complex) matrices, A, satisfying A T R n +R n A = 0. 
These matrices are also said to be perskewsymmetric. Per symmetric matrices are those 
matrices, X, which satisfy X T R n = R n X. P(n,R) (respectively P(n,C)) is the set of 
matrices (real/complex), X, which satisfy X T R n X = R n . 



On In 

• J211 is the 2n x 2n matrix which, in block form, is given by Ji n = 

sp(2n, R) and sp{2n 1 C) represent those real (resp. complex) 2n x 2n matrices which 

satisfy X T J 2n + JinX = 0. Such matrices are also called Hamiltonian. Matrices, Z, 

satsifying Z T J2 n = JinZ are called skew- Hamiltonian. 
( \ 



. so(p, q, R) and so(p, q, C) represent the real (resp. complex) nxn 
matrices (n = p + q), X, satisfying X T I Piq + I p , g X = 0. 



• The anti-trace of an n x n matrix is the sum of the elements on its anti-diagonal. X, 
n x n, is anti-scalar if X = jR n , with 7 e R (or C). 

• Throughout H will be denote the quaternions, while P stands for the purely imaginary 
quaternions, tacitly identified with R 3 . 

Remark 2.1 • i) Throughout this note, use of the following observation will be made: 
Let X be an n x n matrix satsifying X 2 + c 2 I n = 0, c ^ 0. Then e x = cos(c)/„ + 2H^X. 
Here c 2 is allowed to be complex, and c is then taken to be ^/re l i , with c 2 = re 10 ,9 e 
[0,2tt). 

• Occassionally the fact that any matrix which satisfies X 3 = —c 2 X, c ^ 0, satisfies 
e x = I + £HlMx + i-cos(c) x2 Rodrigues - s formu i a ) wU1 also be used . 



• Hi) Explicit formulae for e A can be produced if the minimal polynomial of A is known 
and it is low in degree (cf., [6] where such formulae are written down from the charac- 
teristic polynomial) . Since it is possible to find the minimal polynomial of many of the 
matrices considered here explicitly (i.e., without any spectral information), this removes 
the need for the spectral decomposition, mentioned in the introduction, for A symmet- 
ric. However, since the corresponding explicit formulae for e A are more complicated 
than the ones in i) and ii), they will not be pursued here. See the conclusions for an 
illustration of this issue. 

H®H and gl(4,R): • Associate to each product tensor p®q £ H ®H, the matrix, M Pj9 , of 
the map which sends x e H to pxq, identifying R 4 with H via this basis {1, i, j, k}, Extending 
this to the full tensor product by linearity, yields an algebra isomorphism This connection, 
which is known from the theory of Clifford Algebras, has been put to great practical use in 
solving eigenvalue problems for structured matrices by Mackey et al., [2, 3, 4]. It can also be 
used for finding exponentials, e A , A e gl(A,R) via the following procedure: 

General Algorithm for e A Using H ® H 

• I) Identify u £ H (g> H, corresponding to A via this isomorphism. 

• II) Find e u € H <g> H (in general, this will be possible in closed form only if u (and, 
hence A) possesses additional structure). 

• III) Find the matrix M corresponding to e u - this is e A . 

Note: Throughout this work, tacit use of H <g> H representations of matrices in gl{&, R) will 
be made. These can be easily obtained from the entries of the 4x4 matrix in question (see 
[2, 3, 4] for some instances). In particular, i? 4 = Mj®i, J4 — Mi®j. 

3 Exponentials of Structured 4x4 Matrices 

In this section the algebra isomorphism between H <g> H and real 4x4 matrices will be 
used to find exponentials of various structured matrices. For many of these matrices, their 
exponentials can be found directly from their H®H representations. These will be presented 



first. For the remaining the singular value factorization of matrices, no bigger than 3 x 3, is 
needed. This can be done in closed form, [8]. These will be presented next. 

3.1 Exponentials Directly From H ® H Representation 

Below a (by no means exhaustive) list of nine families of real 4x4 matrices, whose exponentials 
can be directly found from their H <g> H representations, is presented. These families seem to 
be ubiquitous in applications. 

1. 4x4 skew-symmetric matrices: The corresponding element m H®H is p®l + l®q 
with p,q e P. For finding its exponential, it is noted that p <£> 1 and 1 ® q commute, so 
the exponential of the sum is the product of the individual exponentials. Now consider 
Consider (p®l) 2 = — (|| p || 2 )1<X>1. Thus, p®l is annihilated by a quadratic polynomial, 
and its exponential is therefore is [cos(|| p ||)1 + ^T&m^ ff] ® 1 = a; <g> 1, with x a unit 
quaternion. Likewise, e^ 1 ® 9 - 1 = 1 ® y, y a unit quaternion. Thus, e A is the matrix M x ^ y 
- which is a different point of view on the SU(2) x 5/7(2), SO(4, R) relation. 

2.4x4 perskewsymmetric matrices: Such matrices P have H (g> H representations 
p(g>i+a(j(g>l)+j(g>q+f3(l(g>i) = X+Y+Z+W withp G span {i, k},qe span {j, k},a,/3 £ 
R, we find X, Y both commute with each of Z and W. Hence e p = e^ x+y ' e^ z+w \ 
Further, XY = -YX, ZW = -WZ. 

Next, since (X + Y) 2 = (|| p \\ 2 +a 2 )l <g> 1, e {x+Y ^ = cosh(A 1 )(l ® 1) + S ^±±(X + 
Y) = cosh(Ai)(l ® 1) + sin \[ Xl) {p®i + a(j ® 1)), with Ai = y/f}\ p || 2 +a 2 ). Likewise, 
e (z+W) = cosh ( A2 ) (1 i) + EHihl (Z + W) = cosh(A 2 )(l ® 1) + SiH^l (j ® g + /?(1 ® i)), 
with A 2 = VGI g || 2 +/3 2 ). 

Hence, e p is the matrix representation of {cosh(Ai)(l ® 1) + sm ^ Al ) [(^ ® z) + a(j ® 
l)]}{cosh(A 2 )(l ® 1) + 2i2^1 [(j ® g + 0(1 ® i)]}. 

3.4x4 skew-Hamiltonian Matrices: Such matrices, 5, have H ® H representations 
of the form 6(1 ® 1) + p ® j + 1 ® (ci + dk), with b,c,d e R and p £ P. Clearly the 
6(1 ® 1) component commutes with the remaining summands. Thus e s = e b exp(p® j + 
1 ® {ci + dk)). Now, (p® j + 1 ® (d + dk)) 2 = -(- || p || 2 +c 2 + rf 2 )(l <x> 1). Indeed the 
two summands summands anti-commute, while {p® j) 2 =|| p \\ 2 (1 ® 1); (ci + dk)) 2 = 



-(c 2 + d 2 )(l <g> 1). Hence e s = e 6 (cos(A)(l <g> 1) + ^^-(p ® j + 1 <S> (ci + dfc)))> with 
A = ||p || 2 +c 2 + d 2 ). Note A e C. 

4. Five Jordan Algebras: See Appendix II 

5. Eight Lie Algebras: See Appendix II. In particular, one member of this list is precisely 
so(2,2,R). 

6. Simultaneously Hamiltonian, Symmetric, Persymmetric Matrices: These have 
H(g>H representations of the form M = X + Y+Z = f3(j®i)+j(i®k) + 8(k®k), (3, J, 6 e 
R. Now X commutes with both Y, Z while Y, Z anti-commute, and each of X, Y, Z 
squares to a positive constant times 1 ® 1. Hence e M is the matrix representation of 

e x e (Y+z) = [ cos h(/3)(l0l)+ s inh(/3)(i®i)][cosh(A)(l®l)+^^(7(i®fc)+5(fc®fc)],A= ^ffTP 

A 

7. Some Symmetric Toeplitz Matrices: The general case of a symmetric, Toeplitz 
matrix is subsumed by the case of bisymmetric matrices - see Remark (3.3) below. 
Here we identify two important classes which do not require the intervening spectral 
factorization calculations for the general case. 

• Symmetric, Toeplitz, Tridiagonal Matrix: Since such a matrix is met frequently 
in applications, it worth noting that its exponential can be directly computed in 
closed form. Indeed, their H <g> H representations are given by a(l <g> 1) + |(j ® 
i) + \{i ® j) + b(k <g> j),a, b e R. Expressing this as X + Y + Z + W, we see X 
and Y commute with both Z, W and further XY = YX, ZW = -WZ. Hence 

e (X+Y+Z+W) = e X e Y e (Z+W) = e a[ cosh ( g, X ) + S inh(|)(j ® *)] [ C Osh(c)(l ® 1) + 

SS|^((i + A ! )®j)], c =^6. 

• Symmetric, Toeplitz Matrix S Satisfying S13 = 0: This implies that the second 
superdiagonal and subdiagonal vanish. Such matrices have if ® H representations 
of the form a(l <£> 1) + 6(j <8> «) + c(i ® + (X> j). Now, the first and second 
summand commute amongst themselves and with the remaining summands. While 
the third and the fourth anti-commute. Hence, 

e s = e a [cosh(&)(l®l)+sinh(k)(j(^)]^ Vb 2 + c? 



8. Certain Normal Matrices: The general case of normal matrices is subsumed by the 
algorithm below for a symmetric matrix, since the case of skew-symmetric matrices has 
already been dealt with (a matrix is normal iff its symmetric and skew-symmetric parts 
commute). Here we discuss a subclass which does not require the spectral factorization 
calculations needed for exponentiating a symmetric matrix. This subclass is described 
via the following: 

Definition 3.1 Consider a normal TV = S + T, with S its symmetric part and T its 
skew-symmetric part. Expressing T as the sum of two commuting skew-symmetric 
matrices, T\ — M s ®\ and T2 = M\^t,s,t e P it is assumed that || s ||^|| t ||. Such 
matrices will be called special normal. 

Note that special normality forces T^O. Special normality also implies that [S, Tj] = 
0, i — 1,2 (this will be shown below). It is this condition that makes exponentiation in 
closed form possible. 

Indeed, consider, first the case that 7\ ^ 0. Letting S = a(l (g) 1) + p ® i + q (g) j + r ® k, 
the assumption [S 1 , Ti] = forces, in conjunction with the linear independence of the 
elements e x <g> e y , e x ,e y — k, each of the p <£> i,q <E> j,r <B> k to commute with s ® 1. 
This implies that each of p, q, r is parallel to s and hence the symmetric part of iV can 
be expressed succinctly as 

a(l ® 1) + s ® t, s,ie P 
Now the condition, [S, T 2 ] = forces t to be parallel to t. Hence we find 

e N = [e a (cosh(A)/ 2 + (s £)](e s 1)(1 e*), A =|1 « I ' 

A 

If Ti = 0, then the condition [S, T 2 ] = implies that each of p, q 7 r are parallel to one 
another (w.l.o.g p ^ 0), and hence S = p®i, with t = kt,k e R. Hence the above 
formula holds with minor modification. 

Next, it will be shown that special normality implies [S, Tj] = 0, i = 1, 2. One first shows 

T 4 + 2(|| S || 2 + ||t|| 2 )r 2 + [(|| S || 2 -||t|| 2 ) 2 ]/ = (3.1) 

The calculation leading to the above simultaneously shows i) T's minimal polynomial 
is quadratic iff either of s or t vanishes (in this case, trivially [S 1 , Tj] = 0, z = 1, 2; ii) T's 



minimal polynomial is cubic iff || s \\ — \\ t ||. Hence, w.l.o.g T's minimal polynomial is 
quartic, i.e., T is non- derogatory. 

Next, since S, T commute, they are simultaneously diagonalizable, via some unitary 
matrix U. Consider U*TU = U*T{U + U*T 2 U. The last two matrices commute (since 
Ti,T2 commute) and their sum is diagonal. If the entries of the diagonal matrix U*TU 
are all distinct and non-zero, then the matrices U*T\U and U*T2U are themselves 
diagonal. Thus they also commute with U*SU, which implies that S commutes with 
both the Tj. Note T being non-derogatory implies the assumptions about the diagonal 
entries of U*TU, in view of the nature of the eigenvalues of a 4 x 4 skew-symmetric 
matrix. 

9. Certain Non-Toeplitz Bisymmetric Matrices Every persymmetric matrix is of the 
form RS, with S symmetric. Similarly, Hamiltonian matrices are of the form JS, with 
S symmetric. Such matrices can often be exponentiated in closed form, if in addition, 
R±S = SRa, (resp. J 4 5 = SJ A ). 

Indeed, since RS = SR, and R 2 = I, we find e RS = cosh(S) + i?sinh(S) (this equation 
is valid in any dimension). Now, S = a(l <g> 1) + p® i + g® j + r ® k satisfies R4S — SR4 
iff i) p is parallel to j and ii) q,r are perpendicular to j. If, in addition we suppose 
either q, r are parallel or q, r are perpendicular to one another, then exponentiation in 
closed form is possible. For brevity the former possibility is assumed. Hence 



Note, in particular, that RS is symmetric, persymmetric, but not Toeplitz. 

Writing RS as R(p,h + S), with S = X + Y, we sec that it suffices to find e RS '. 
Now, notice that X and Y commute and X 2 = e 2 I,Y 2 = (a 2 + /3 2 )(j 2 + S 2 )I = \ 2 I. 
Hence, cosh(^) = cosh(X) cosh(F) + sinh(X) sinh(Y), and sinh(5) = sinh(X) cosh(Y") + 
sinh(Y)coshpO. Butsinhpf) = s -i2^il X ; sinh(Y) = SiH^i Y; cosh(X) = cosh(e)/; cosh(r) = 
cosh(A)/. Hence e RS is the matrix given by: 



5 = o(l <g> 1) + e(j ® i) + (ai + /3k) <g> (jj + 8k) 



[cosh(/x)/4+ 



sinh(/i) 
A* 



R] [cosh(e) cosh(A) 1+ 



sinh(e) sinh(A) 



Ac 



XY]+R[ 



sinh(e) cosh(A) 



c 



X+ 



sinh(A) cosh(e) 



A 



Y] 



Similarly, if JS — SJ, one finds (since J 2 = —I) that 

e JS = cos(S) + J 2n sin(S') 

Now if S, symmetric, commutes with J, then fortunately (or unfortunately) J4S is also 
simultaneously skew-symmetric, and therefore the previous formula is yet another way 
of exponentiating J4S. Hence, the details are omitted. 

3.2 The General Symmetric Case 

Exponentiating the general 4x4 symmetric matrix requires the spectral factorization of a 
3x3 matrix (which can be done in closed form). Before getting to that, the principal enabling 
feature of the algorithm below is described by the following: 

Proposition 3.1 The exponential of a(l <X> 1) + Y^i=i u i ® G P, with {ui,i = 

1, . . . , 3}, {vi, i = 1, . . . , 3} each an orthogonal triple in R 3 is given by e a H 3 =1 e( Ui ® Vi \ with 
e (ui9vt) =cosh(|| Ul HI! Vl ll)(l®l) + 8i g|iBf (uiQvi). 

Proof: It suffices to observes that each of the summands in a(l ® 1) + X^i=i u i ® v i 
commutes with each other due to the orthogonality property. The formula for e ^ Ui ^ v ^ is now 
just a consequence of (ut <g> vi) squaring to a positive constant times the identity. 

Remark 3.1 If the triples {ui},{vi} were instead each parallel to each other, then once 
again the exponential of a(l <£> 1) + X)i=i u i ® v i 1S quickly computed, since now once again 
each summand commutes with one another. There are other possible configurations which 
will render the calculation of the exponential in closed form too. However, these will not be 
pursued here for brevity. 

Remark 3.2 • i) Consider the element p<g)i + q<E)j + r<E> k,p, q, r e P. Then, as observed 
in [2], if Y^i=i a i u i v I ' ,Ui,Vi e R 3 is the singular value factorization of the real 3x3 
matrix, \p \ q \ r], it follows that p ® i + q® j + r ® k = Y^=i a i u i ® v ii where the 
vectors Ui 7 Vi have been identified with the corresponding pure quaternions (in lieu of 
the elegant proof in [2] , one can also verify this via direct calculations which show that if 
for pi, qt, n, Si e P, i = 1, . . . , 3, the 3x3 matrices X)i=i PilJ j Si=i r i s I coincide, then 



Si=i M Pi ® qi — Y^=i M ri ® Si )- Since the {ui\, {vi} are each an orthonormal triple, the 
exponential of p <£> i + q <X> j + r <X> fc, which equals the exponential of X)i=i CT i M j ® can 
be explicitly found by using Proposition (3.1). The only issue is computing the singular 
value factorization of a real 3x3 matrix. However, this is the spectral factorization of a 
real 3x3 symmetric matrix, which itself can be done in closed form, [8] . It is interesting 
to note that the technique described in [8], consisting of 3 x 3 matrix manipulations, 
can itself be implemented via quaternions. 

• ii) Note the subsequent rotations employed in [2] to diagonalize a symmetric matrix are 
not required, since diagonalization is not being employed here to compute exponentials. 
Only the reduction to form used in Proposition (3.1) is needed. 

This leads to the following algorithm for the exponential of a 4 x 4 symmetric matrix: 

• Represent the matrix as a{l ® I) + p ® i + q® j + r ® k,p, q, re P. 

• Compute the singular value factorization, Y^i=i a i u i v J \ u ii v i <= R 3 OI tnc rca l 3x3 
matrix, \p \ q \ r]. 

• Compute the exponential of a(l <£> 1) + Y^i=i a i u i ® v i vl& Proposition (3.1). The 4x4 
matrix representing this element of H ® H is e A . 

Remark 3.3 The special classes of 4x4 bisymmetric matrices (i.e., simultaneously symmetric 
and persymmetric) and 4x4 symmetric and Hamiltonian matrices are, of course, subsumed 
by the foregoing algorithm. However, it is worth pointing out, in view of their importance 
in applications, that the singular value factorization needed is easier to find than in the fully 
symmetric case. Indeed, a bisymmetric matrix is represented by a(l ® 1) + b(j ®i)+p®j + 
q ® k,p,q e span {i, k}, a, 6 € R. Thus, it suffices to find the singular value factorization of 
the 2x2 matrix \p | q] - which is the spectral factorization of a 2 x 2 real symmetric matrix. 
Likewise, a symmetric, Hamiltonian matrix is represented by g <£> i + r <x> k. Thus, it suffices 
to find the singular value factorization of the 3x2 matrix [p \ q] (only two of the left singular 
vectors are needed) . There are many other cases of symmetric matrices possessing additional 
symmetry which are susceptible to the same observation. 



Remark 3.4 Extension to Complex Matrices: Some of the procedures extend to special 
classes of complex matrices. This is illustrated for matrices in so(4, C). Such a matrix can be 
represented in the form a 1 M im +f3 1 M ]m +^iM m i+a2M 1 ®i+(52M 1 ® j +^ 2 M 1 ® k = X)j=i X U 
with ai,fti,'fi G C. Now the fact that these constants are complex does not prevent each 
of Xi,. . . ,X 3 from commuting with each of X4, . . . ,Xq. Neither does it prevent each of 
X\, . . . , X3 anti-commuting with one another nor each of X4, . . . , Xq anti-commuting with 
one another. Finally, Xf = —c-I^ for each i — 1, ...,6, for some e C. Hence the 
exponential is given by 

[cos(Ai)/ 4 + bm f Al \ aiM^ 1 +/3iM J -® 1 +7i^ 

Ai A2 

with Xf = — (a? + 01 + 7 4 2 ),i = 1,2. Similarly the technique for p(4 7 R) extends verbatim 
to p(4, C). However, while the methods based on the singular value factorization extend 
verbatim for purely imaginary symmetric matrices, they are not applicable to general complex 
symmetric matrices. To see what is needed for the extension, consider traceless symmetric 
matrices (w.l.o.g). Let Ar and Ai be the real and imaginary parts of A. Since these are 
symmetric as well, one can associate two triples (Pi,Qi, rj) G P 3 , i = 1,2. Let Mj = [pi \ | r»] 
be the corresponding real 3x3 matrices. If these could be simultaneously brought into the 
canonical forms Mi — J2k=i a k u k v ki with the Uk and Vk orthonormal, a\ E R, then clearly 
the algorithm for real symmetric matrices would extend verbatim to such matrices. Many 
sufficient conditions are known for such simultaneous canonical form, [9] . One such condition 
is that both MiM% , Mi should be symmetric. 

4 Conclusions 

In this note, closed form formulae are provided for exponentials of several important families 
of real (and complex) 4x4 matrices. In conjunction, with techniques such as Givens rotations, 
these formulae provide algorithms for exponentiating classes of structured matrices in higher 
dimensions. The principal technique is the invocation of the associative algebra isomorphism 
between gl (4, R) and H ® H. It is the ease of multiplication in H £g> H which facilitates the 
discovery of closed form exponentials for many matrices 

It is possible to write down exponentials of matrices once their minimal polynomial is 



known (especially if they are at most quartic). However, these formulae themselves can be 
quite complicated and hence they were not pursued in this note. This is exemplified by generic 
4x4 skew-symmetric matrices, whose minimal polynomial is quartic. The corresponding ex- 
ponential formula, though equivalent to the one given here, is substantially more complicated. 
In our opinion most 4x4 matrix calculations should be done in H ® H . The formulae for 
the minimal polynomial of a 4 x 4 skew-symmetric matrix [see Equation (3.1)], without any 
spectral information, is yet another vivid illustration. 

Clearly, H®H is associative algebra isomorphic to <?Z(4, c), where H is the complcxification 
of H. One can identify the latter with gl(2,C). However, it is better to view its elements 
as q = x + x\i + x 2 j + x 3 k,Xi e C and define q = x~o — x~\i — x 2 j ~ x~sk. This notion 
of conjugation is equivalent to Hermitian conjugation in gl(4,C). This does not, however, 
render calculating exponentials in su(4) as simple as in so(4, R) (after all one cannot run 
away from the curse of dimensionality by such an artifice). However, several Hermitian and 
skew- Hermitian matrices (e.g., whose real and imaginary parts come from special normal real 
matrices) are easily exponentiated. 

5 Appendix I: 

In this appendix, a different approach to the exponentiation of matrices inp(4, R),so(2, 2, R),p(3, R) 
is described, which reduces the problem to the exponentiation of 2 x 2 matrices (this is equally 
applicable to their complex counterparts). This is first illustrated for matrices in so(3, R) and 
so(4, R) since this should be reasonably well known terrain. Attention, in particular, is drawn 
to Remark (5.1), which provides the correct heurisitics needed to generalize this to the ma- 
trices in p(4, R), so(2, 2, R),p(3, R). 

Consider an element A e so(3,R). Its exponential can be computed explicitly via the 
Rodrigues formula. The usual derivation of this relies on the fact that A satisfies 

A 3 + X 2 A = 0, A e R 

Any matrix which satisfies this equation will satisfy the Rodrigues formula. There is a 
equally well-known relation between su{2) and so(3,R). What is, perhaps, less appreciated 
is that this relation yields an explicit technique to find e A ,A e so(3,R). To describe this, 



fix a G £ SU{2). Consider V = {A \ A* = A,Tr(A) = 0}. SU{2) acts via conjugation on 
elements A £ V, viz., 4>g{A) — GAG -1 . It is well known, that upon identifying V with R 3 
through the basis {<Jk,k — x,y,z}, this action yields a proper rotation of R 3 . Thus, we get 
a homomorphism, : SU{2) — ► SO(3,R), which sends G to the matrix of <pa with respect 
to the basis {<7fc, k = x,y, z}. This is a surjective, two-one, homomorphism. Linearizing this 
map, we get a Lie-algebra isomorphism ip : su(2) — > so(3, J?), viz., ^>(A) is the matrix of the 
linear map which sends w £ V to Aw — vA with respect to the {a k , k — x,y, z} basis, with 
A £ su(2). This is a Lie-algebra isomorphism. From elementary considerations in Lie theory 
-0 and 4> provide the following technique to find e A ,A£ so(3, R): 

• i) Find B = ^(A) in su(2) 

• ii) Compute e B £ SU(2) - this can be explicitly done since satisfies the condition in i) 
of Remark (2.1)., 

• iii) Compute the matrix <p e B - this is e A . 

This is arguably easier to use than the Rodrigues formula (it is left to the reader to verify 
that the two result in the same formula) . This is not to disparage the Rodrigues formula - it 
applies to situations where Lie theory would have no visible role. But the fact that a 3 x 3 
exponential has been computed with a 2 x 2 calculation is significant. Similar and even better 
savings occur by such arguments. 

Exponentials in so(4, R): There is a well known two-one Lie group homomorphism denoted 
by (f> : SU{2) x SU{2) -> 50(4, R), given by the action of SU(2) x SU(2) on the vector space, 
V, of real linear combinations of I 2 ,iai-,k = x,y,z, viz., for fixed G,H £ SU(2) x SU{2), 
let 4>g,hV — > V be given by (Pg,h(X) — GXH^ 1 ,X £ V. Once again this is a proper 
rotation of R 4 (identified with V via this basis), and <j>(G,H) is precisely the matrix of 
this map with respect to this basis. Linearizing this gives a Lie algebra isomorphism, -0 : 
su(2) x su(2) — > so(4, R) which sends (X,Y) £ su(2) x su{2) to the matrix of the map (with 
respect to the Ii,ia k basis) which sends Z £V X,o XZ — YZ. This yields an algorithm to find 
e A , A £ so(4:,R), which reduces to finding two 2x2 exponentials in su(2) - the statement of 
the algorithm is omitted (mimick the p(4, R) algorithm given below). 



The corresponding relations between SL(2, C) (respectively SL(2, C) x SL(2, C)) and 
SO(3, C) (respectively 50(4,(7)) once again reduce exponentiation of matrices in so(3,C) 
and so(4, C) to 2 x 2 calculations. Note that the fact that 50(3, C) etc., are not compact 
docs not matter for the veracity of this procedure. All that is needed for finding e A is that 
the corresponding be a Lie group homomorphism (it need not even be surjective) and the 
corresponding -0 be a Lie algebra isomorphism. 

Remark 5.1 Traditional proofs of the SU(2) covering of SO(3,R) proceed by i) using 
su(2) itself as the vector space V, and ii)then, by constructing a bilinear form, K(X, Y) = 
Tr(ad X ad Y) on su(2) and showing that this is preserved by the action of SU{2). For our 
purposes it is more useful to proceed differently. On any (sub)space of 2 x 2 matrices, there 
are two obvious candidates for quadratic forms, viz., i) Tr(X 2 )- and ii) det(X). One is even 
lead inexorably to these forms upon inspecting the forms of the maps <j) used above for both 
so(3, R) and so(4, R). Polarizing these two leads to the following choices: 



It is easy to see that, with the choice of bases made in the derivation of the so(3, R) 
(resp. so(4, R)) algorithms, that the symmetric matrices representing these two forms are 
precisely 2I 3 (resp. I4). This immediately shows that the matrix of the corresponding 0's are 
orthogonal. 

Remark 5.2 Lorenz Lie Algebra: Here a different perspective on the work of [7] on the 
exponentials of matrices in so(l,3, R) is provided. Indeed, letting V be the i?-linear span 
of {I 2l (Tie, c y , <r z } (i.e., V is the space of 2 x 2 Hermitian matrices), it is found that the 
matrix of L 2 (X, Y) is precisely 27i i3 . If SL(2, C) acts on V via <j> M {v) = MvM*,v e V, M e 
<SX(2,0), then L 2 (X,Y) is preserved and the matrix of <Pm in this basis is in the Lorenz 
group. Linearizing <f), we get a technique to find exponentials in so(l,3,R), cf., [7]. 

Below the same thinking is used to compute exponentials in p(4, R), so(2, 2, R) and p(3, R). 
The method can be applied to several other Lie algebras stemming from symmetric, non- 



L 1 (X,Y) = Tr(XY) 



(5.2) 



L 2 (X, Y) = det(X + Y) - det(X) - det(Y) 



(5.3) 



degenerate, bilinear forms on R 4 . However, we limit ourselves to these cases for brevity. 

Exponentials in p(4, R): 

Consider gl(2, R), identified with i? 4 via the basis, {En, E\ 2 , —E 2 i,E 22 }. Let SL(2, R) x 
SL(2,R) act on gl(2,R), via 4>g,h(X) = GXH^ 1 . This action leaves the bilinear form 
L 2 (X, Y) of Equation (5.3) invariant. Furthermore the symmetric matrix representing it, with 
respect to this basis, is precisely i?4- Thus the matrix of <Pgm is in P(4, R). Linearizing this 
we get a Lie-algebra isomorphism (that this is a Lie- algebra homomorphism is standard - it is 
easily verified that it is an isomorphism): ip : sl(2, R) x sl(2, R) — + p(4, R), which sends a pair 
(g, h) £ sl(2, R) x sl(2, R) to the matrix of the linear map L 3yh (X) = gX - Xh, X e gl(2, R) 
with respect to the {En, E\ 2 , —E 2 i,E 22 } basis. This leads to the following algorithm to find 
e A ,4ep(4,4 

Algorithm for e A ,Ae p(4, R): 

• i) Find the pair (g,h) = ^(A) e sl(2,R) x sl(2,R); 

• ii) Find G = e 9 , H = e h . This is easily done since g, h satisfy the equation in Remark 
(2.1 -i)) 

• hi) Find the matrix of <Pgm with respect to the above basis. This is e A . 

Exponentials in so(2,2,R): 

Now identify gl(2,R) with R 4 via the basis {I 2 ,Ei 2 — E 21l cr Xl a z }. Then the matrix of 
L 2 (X,Y), with respect to this basis is precisely I 2 , 2 . Let SL(2,R) x SL(2,R) act on gl(2,R) 
via <Pgm(V) = GVH-\G,H e SL(2,R),V E gl(2,R). Then (j> GM preserves L 2 (X,Y) and 
hence its matrix, with respect to this basis, is in SO(2,2,R). Linearizing this action we get 
a Lie algebra isomorphism tp : sl(2, R) x sl(2, R) — > so(2, 2, i?), with tp(g, h) being the matrix 
of the linear map ip g ,h(v) = gv — vh, v e gl(2, R) with respect to the same basis. This leads 
to an algorithm, similar to the previous one, for finding e A , A e so(2, 2, R). 

Exponentials in p(3,R): 

Now identify R 3 with the real span of E\ 2l o x ,E 2 \. This is sl(2,R). Then the matrix 
of Li(X,Y), with respect to this basis, is, upto a constant, R$. Let SL(2,R) act on this 
copy of R 3 via (f>G{h) — GhG^ 1 . This action preserves L\(X,Y). Thus, the matrix of 



(j)c is in P(3, R) and the map (f> : SL(2, R) — » P(3, i?) is easily seen to be a Lie group 
homomorphism. Linearizing <j> leads to a Lie algebra isomorphism ip : sl(2,R) — > p(3, i?) 
which sends ft G sl(2,R) to the matrix of the linear map, which sends X e sl(2,R) to 
/iX — X/i (identifying sl(2,R) with i? 3 via the above basis). This leads to an algorithm for 
finding e A ,Ae p(3,R). 

Remark 5.3 i) The last calculation can be mimicked to find exponentials in so(2, 1,R). 
Indeed, identify sl(2, R) with R 3 via the basis {a Xl cr z , E\ 2 — E 2 {\ and proceed verbatim as in 
the p(3, R) case, ii) All of the above calculations extend to find exponentials in p(4, C) etc., 
The only difference is one works with complcxifications of the various Lie algebras introduced 
before, i.e., sl(2,C) x sl(2,C) for p(4,C) etc., 

6 Appendix II 

In this appendix are listed i) five classes of matrices, each a Jordan algebra, which can be 
exponentiated by mimicking the technique for skew-Hamiltonian matrices; ii) eight classes of 
matrices, each forming a Lie algebra, which can be exponentiated by mimicking the technique 
for perskewsymmetric matrices. In most cases the technique extends to their complex ana- 
logues (e.g., so(2,2,C)), cf., Remark (3.4). In both lists, both the H ®H representation and 
the 2x2 block representations are provided. 

Remark 6.1 Let M\,M 2 be two invertible, symmetric (resp. skew-symmetric) matrices, 
with the corresponding bilinear form on R n denoted by <, >m x , <, >m 2 ■ The two forms are 
defined to be equivalent if there is an orthogonal matrix G such that G T M^G = M 2 . If this is 
the case then the corresponding Jordan algebras, Jj = {X \ X T M,- L = MiX}, i = 1, 2 and the 
corresponding Lie algebras Lj = {X \ X 7 'Mi = —MiX},i = 1,2 are conjugate. Specifically 
J 2 = G T JiG,L 2 = G T L\G. Thus, if one knows exponentials of matrices in Ji (resp. L\), 
then one can find exponentials of matrices in J 2 (resp. L 2 ) provided G is explicitly described. 

In the first list, the first two Jordan algebras pertain to bilinear forms which are equivalent 
to J 4 , while all the matrices in the second list stem from symmetric forms equivalent to i? 4 . 
While it is possible to explicitly construct the corresponding G's, it is far easier to work with 
the matrices in these lists directly. 



Exponentials of Five Jordan Algebras 



p® A; + a(l(g)l) + l® (bi + cj), a, b, c € R,p ^ P. The block representation is 
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i(E)<7 + a(l(g)l) + (6j + cfc)(Xil,a,6,ce R 7 q& P. The block representation is 



with A, D some 2x2 scalar matrices, B = 
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k <g> q + a(l ® 1) + (bi + cj) ® l,a, b, c G R, q G P. The block representation is 
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Exponentials of Eight Lie Algebras 



V 7 6 J 



• so(2, 2, R). The H <8> H representation is a(l <£> i) + p <g> i + b(i ® i) + i <X> q,p, g € 

' .4 ^ 



span {j, fc}, a, 6 e i?. The block representation is 



B C 



, where B is any 2x2 



matrix, while A, C are 2x2 anti-diagonal matrices with zero anti-trace. 



p <8> j + a(j ® 1) + j ® q + b(l ® j),p.g € span {z, fc}, a, 6 € -R. The block representation 
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where B is any 2x2 matrix, while A, C arc 2x2 anti-scalar matrices. 
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p <g> i + a(k <g> 1) + k <g> g + 6(1 <g> i),p G span {i, fc}, g G span {j, fc}, a, 6 G i?. The 
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p <8> j + a(k ® 1) + k <g> g + 6(1 <g> j),p € span {i, j}, g G span {i, fc}, a, 6 G i?. The block 

, with A, D zero-trace, diagonal 2x2 matrices, while 



representation is 
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p ® j + b(i <g> 1) + a(l <g> j) + i <g> g,p G span {j, fc}, g G span {i, fc}, a, 6 G i?. The block 



representation is 
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p®k + a(i®l) + b(l®k)+i®q 7 p£ span {j, k}, q G span {i, j}, a, 6 G _R. The 2x2 block 



representation is 
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, with B,C 2x2 zero anti-trace, anti-diagonal matrices, 
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p <8> k + a{j <X> 1) + b(l <X> k) + j <X> q,p € span {i, fc}, g G span {i, j}, a, 6 G _R. The block 

, with A, D 2x2 diagonal, zero-trace matrices, while 



representation is 
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